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Abstract: Due to instanton effects, gauge-theoretic large N expansions yield asymptotic series, 
in powers of 1/A^^. The present work shows how to generically make such expansions meaningful 
via their completion into resurgent transseries, encoding both perturbative and nonperturbative 
data. Large N resurgent transseries compute gauge-theoretic finite N results nonperturbatively 
(no matter how small N is). Explicit calculations are carried out within the gauge theory 
prototypical example of the quartic matrix model. Due to integrability in the matrix model, it 
is possible to analytically compute (fixed integer) finite N results. At the same time, the large 
N resurgent transseries for the free energy of this model was recently constructed. Together, 
it is shown how the resummation of the large N resurgent transseries matches the analytical 
finite N results up to remarkable numerical accuracy. Due to lack of Borel summability, Stokes 
phenomena has to be carefully taken into account, implying that instantons play a dominant role 
in describing the finite N physics. The final resurgence results can be analytically continued, 
defining gauge theory for any complex value of N. 
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1 Introduction and Discussion 

The large N expansion of nonabelian gauge theories was hrst uncovered four decades ago [1] 
(see, e.g., [2] for a modern review). On the one hand, this expansion is topologically organized 
as a genus expansion which has led to remarkable developments in large N duality, e.g., from 
supersymmetric gauge theories [3] to broad classes of matrix models [4]. On the other hand, 
this expansion may be regarded as an approximation to compute gauge-theoretic observables 
nonperturbatively (from the point-of-view of the gauge coupling) at low values of the rank of the 
gauge group^. While the first spin-off has been extremely successful over the years, the second 
has not received much attention beyond the planar approximation and its first few (mainly 
perturbative) corrections. It is our goal in this paper to introduce and develop a rather general 
method to extract (eventually exact) finite N results out of the large N expansion. 

The reason why at first sight the large N expansion seems meaningless is because the resulting 
genus expansion is an asymptotic series, i.e., it has zero radius of convergence. Consider the free 
energy F = log Z of some nonabelian gauge theory, and denote by t the’t Hooft coupling (held 
fixed in the ’t Hooft limit). If one starts at fixed genus, say with a perturbative calculation of 
the planar limit Fo(t), this results in a convergent series in t with non-zero radius of convergence 
[5, 6]. But when one starts addressing the genus expansion of the free energy, one finds that 
its genus g perturbative coefficients, Fg{t), grow factorially fast as Fg ~ {2g)\. This results in a 
divergent asymptotic series which is not even Borel summable [7, 8] . How can one make sense of 
such an expansion, in order to obtain finite N results out of the large N expansion? 

For concreteness let us focus upon the free energy, whose large N asymptotic expansion is 

+ 00 

( 1 - 1 ) 

5=0 

One might hope to make sense out of this factorially divergent series via Borel resummation. In 
fact, the Borel transform removes the factorial growth from the perturbative coefficients as 

„2a—3 

^As originally emphasized by’t Hooft [1], it was then hoped that the large N expansion, and in particular the 
planar limit, would lead to “reasonable” approximations in spite of the physical interest in smaller values of N. 
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so that its (non-zero) radius of convergence is now dictated by the subleading exponential growth 
of the original perturbative coefficients. Upon analytic continuation of i3[U](s) along the complex 
Borel s-plane, the Borel resummation follows whenever the ray of integration along the direction 
6 avoids the integrand’s singularities. In this case one finds 

^e‘®oo 

S9F{N,t)= dsB[F]is)e-^^. (1.3) 

Jo 

The important point is that even when the choice of direction 6 is such that the above integration 
is possible, i.e., that one finds Borel summability along 6, this expression need not match the 
exact result due to the occurrence of Stokes phenomena [9]. In fact, the factorial large-order 
behavior was already telling us that one is missing nonperturbative corrections, typically of the 
form ~ exp(—A^), which are invisible to the perturbative expansion (1.1). But in spite of that, 
these initially exponentially suppressed terms may not be neglected. In fact, upon variation of 
9 in (1.3), they may grow to become of order one or, eventually, exponentially enhanced with 
respect to (1.1), and thus completely obliterate any sense in which the resummation (1.3) may 
yield correct results. This is the essence of Stokes phenomena, implying that the perturbative 
series alone cannot properly define the free energy and it needs to be enlarged into a transseries. 

That Borel resummation, if allowed, cannot be the whole story when it comes to the genus 
expansion was recently verified numerically to great precision in [10] (following earlier work on 
the quantum mechanical quartic oscillator [11, 12]). In the example of the large N ABJM gauge 
theory (but having also addressed other examples, such as topological strings in the local x 
Calabi-Yau geometry), [10] showed numerically that even if this case satisfies the requirements of 
Borel summability, Borel resummation of the genus expansion does not agree with the exact value 
(computed via integrability), with the mismatch controlled by complex D2-brane instantons. 
As explained in the above paragraph, this is just explicitly verifying that, in general, Borel 
resummation is somewhat useless from a practical point-of-view due to Stokes phenomena: even 
if there are no obstructions to performing the inverse Borel transform in (1.3), and the asymptotic 
expansion is dubbed Borel summable, the resulting resummation will yield the correct result 
only if Stokes phenomena has not yet occurred. Otherwise it will always miss nonperturbative 
contributions: further terms in the transseries may have already been turned on by the crossing 
of a Stokes lines and these are thus needed and crucial to obtain correct results. 

Thus, building upon work initiated in [13] and as we shall discuss in the main body of the 
paper, finite N results can only be obtained out of transseries; schematically of the form 


+O0 H-OO ^ +O0 - 

F{N, t) = J2 (t) + ^ ^ fW (t) + ^2 ^-2NA(t) ± Ff) {t) + • 

g=0 g=0 g=0 


(1.4) 


Here A{t) is the large N instanton action^ and Fg^\t) the perturbative ^-loop amplitudes around 
the nonperturbative n-instanton sector. Generically all these perturbative series are asymp¬ 
totic. At large N, the exponential terms are non-analytic (thus the term “transseries”), and the 
transseries parameter a encodes Stokes phenomena (we shall discuss this point carefully in the 
main body of the text, but see also [18] for general and very explicit formulae). Perhaps one of 
the most important properties of these transseries is that they are resurgent. This means that, 
chosen any perturbative or multi-instanton sector, the coefficients of its associated asymptotic 

^In those particular cases where the gauge theory is given by a matrix model, then the instanton action has a 
very physical description given in terms of eigenvalue tunneling; see [14-17]. 
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series, Fg'^\t), have their large-order growth precisely dictated by the instanton action and by 
(the low-orders of) “nearby” coefficients F^\t). Conversely, decoding the asymptotics of a 

chosen sequence Fg^\t) into its different power-series (analytic) and exponential (non-analytic) 
contributions in one will find the “resurgence” of all other sectors with n' n, g' g. See, e.g., 
[19] for very explicit large-order resurgence formulae. All distinct components of the transseries 
(1.4) are thus very tightly bound together by a stringent web of resurgence relations. 

Some light introductions to resurgent transseries may be found in, e.g., [20-22] with more 
mathematical flavor, or in the first sections of [19] with a flavor closer to the present work. 
Resurgent transseries are at the basis of Ecalle’s work [23], where by joining analytic and non- 
analytic building blocks they allow for representations (in some sense, for constructions and 
definitions) of broad classes of functions, including functions with singularities and branch cuts. 
This is why in many cases where one might only have large N expansions available, a large N 
resurgent transseries may allow for a definition of the finite N theory and (considering arbitrarily 
large numbers of components in the transseries) for a calculation of exact finite N results. 

While the initial gauge theory was only defined for integer values of N, the resurgent 
transseries is constructing a function of N, eventually valid upon analytic continuation for arbi¬ 
trary complex values of N. While at first this might sound unexpected for nontrivial interacting 
gauge theories, it is actually familiar in free gauge theories such as the Gaussian or Penner matrix 
models (see, e.g., [24, 25] for brief reviews). For instance, the partition function of the Gaussian 
matrix model (normalized against the volume of its U(A^) gauge group), while only defined for 
integer values of N via a matrix integral, may still be computed as^ 

rF_ 

Zg{N) = ^^G2{N + 1), (1.5) 

(27r)2 

where G 2 {z) is the Barnes G-function which is in fact an entire function on the complex plane 
(see, e.g., [26]). Finding nontrivial gauge-theoretic examples where nonperturbative completions 
allow for results defined at arbitrary values of N is not as straightforward. But a class of very 
interesting such examples was recently addressed in [27], within the context of ABJM gauge 
theory on the three-sphere. In their examples, supersymmetry and integrability ensure that all 
(genus) contributions to the partition functions actually truncate to finite sums, allowing for 
the computation of closed-form expressions at arbitrary values of N. In particular, partition 
functions for these theories were also found to be entire functions on the complex Wplane. 

One fascinating spin-off of the calculations in [27] is the possibility to study exact prop¬ 
erties of the partition function of ABJM gauge theory as a function of A^. As shown therein, 
the fact that all partition functions found in [27] were entire functions on the complex iV-plane 
implies the non-existence of phase transitions in iV; in particular the non-existence of a phase 
transition as one moves from (real) large N to small N. From a dual holographic viewpoint 
[3] this implies there is really no drastic distinction between a “quantum” small-distance and a 
“semiclassical” large-distance phase, at least in this class of theories. Furthermore [27], this is 
particularly striking given that if one just looks at the genus expansion by itself, the perturbative 
free energies do have branch cuts in the ’t Hooft coupling t, misleading us on the understanding 
of the nonperturbative physics. This phenomenon where the perturbative singularity structure 
is nonperturbatively smoothened was also found earlier in [28], within the context of minimal 

®The so-called string coupling Qs appears as the overall coupling in front of the quadratic potential inside the 
matrix integral defining the partition function, and relates to the’t Hooft coupling as t = QsN. 
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string theory. Therein, certain observables such as the FZZT partition function, while having 
semiclassical branched domains, end up having entire analytical properties when addressed non- 
perturbatively. Interestingly enough, the mechanism by which this occurs is again based upon the 
aforementioned Stokes phenomena, whereby (initially) exponentially small contributions grow to 
take dominance and dramatically change the structure of the solutions. The large N resurgent 
transseries methods we develop in this paper will allow for these types of analysis in broad classes 
of gauge theories, even those where neither supersymmetry nor integrability are present. 

One final point to discuss is whether the resurgent transseries completion we propose, to 
non-integer (and complex) values of N, is the only possible such completion? In principle, gauge 
theory observables are only defined at positive integer values of N, and, clearly, many possible 
interpolating functions may exist through such a discrete set of points. It was already suggested 
in [27] that the familiar Gamma function serves as a good analogy. Certainly there are many 
other functions which solve the factorial interpolation problem, and which in fact behave much 
better (the Gamma function is meromorphic). One such example is the Hadamard Gamma 
function, which interpolates the factorials as an entire function, but many others exist; see [29] 
for an excellent review. What tells all these possible interpolating functions apart is, of course, 
whatever (differential or functional) equation defines them. In this regard, the usual Euler 
Gamma function is the only^ function which satisfies the continuous (functional) version of the 
factorial recursion, 

T{z + 1) = zT{z). (1.6) 

When translating to our current large N problem, a rather similar situation still holds. As we 
shall see in the main text of this paper, the partition function at hnite integer N is computed by 
solving a discrete recursion relation (the so-called string equation [30], but we shall discuss this 
at greater length later), while the large N resurgent transseries is obtained out of the continuous 
(functional) version of this discrete string equation [13, 19, 31]. Furthermore, this solution is 
very special, in the sense that it is a resurgent solution: all its perturbative and nonperturbative 
building blocks are very tightly bound together! Although we do not address this point in this 
work—and it should certainly be addressed in future research—it might be reasonable to expect 
that these two properties combined provide a unique analytic continuation of gauge theoretic 
observables to complex N (at least up to some choice of boundary conditions). 

The problem of resummation of divergent asymptotic series in Physics has a long history; 
see, e.g., [32]. In order to apply it, and extend it via the use of resurgent transseries, to the 
problem of the large N expansion one hrst needs data, on both perturbative and multi-instanton 
expansions. To the best of our knowledge, the gauge theories for which at present there is 
the largest amount of available large-resurgent-transseries data are matrix models (and their 
double-scaling limits); in particular matrix models with polynomial potentials [13, 33, 34, 19, 31]. 
This type of matrix models is a very reasonable prototype for generic gauge theories. We shall 
perform our explicit calculations in the example of the quartic matrix model [19, 31], although 
we should stress that our methods are completely general within gauge theory. 

This paper is organized as follows. We begin in section 2 with a calculation of exact results 
in the quartic matrix model, at small (integer) values of N. This is done via standard methods 
using orthogonal polynomials, and several exact results are obtained for the partition function 
and the recursion coefficients. At fixed iV, these quantities are still functions of the coupling, or, 
equivalently, of the’t Hooft coupling, and we address their corresponding non-trivial monodromy 
properties. In section 3 we then address the gauge-theoretic large N resurgent solution to our 

^To be correct, uniqueness is only achieved after adding the technical requirement of logarithmic convexity [29]. 
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example. While the resurgence data of the quartic matrix model is reviewed in appendix A, 
section 3 shows how to make use of this data (or resurgence data from any other gauge theory) 
in order to address the resummation problem and produce hnite N results. Comparison with the 
analytical results is achieved up to remarkable numerical accuracy. Finally, section 4 discusses 
Stokes phenomenon and how the full resurgent transseries is crucial to reproduce the analytical 
finite N results across the complex t-plane, including the distinct non-trivial monodromy prop¬ 
erties one found earlier for different integer values of N. We conclude with a discussion of how 
the large N transseries may be analytically continued for any complex value of N, and what are 
the properties of the resulting partition function on the complex A^-plane. 

We believe that many interesting problems may be addressed in the near future using the 
methods we introduce and develop in this paper. For instance, one may consider gauge theories 
in higher dimensions. When considering supersymmetric gauge theories on compact manifolds 
where some observables localize [35], we may fall back into the realm of matrix models. The 
resurgent analysis of localizable observables in some three and four dimensional gauge theories 
was initialized in [36] and it would certainly be interesting to address hnite N calculations in 
that realm. If we move away from matrix models, the main difference will be that the (now 
diagrammatic) computation of perturbative and multi-instanton coefficients will be much more 
time consuming. Our uses of large N resurgent transseries, however, should hold step by step. 
Perhaps a good route where to start would be to follow the recent resurgence calculations for the 
CP^ model in [37-39] and address the corresponding large N limit in that context. 

Via large N duality, another class of theories to explore are string theories. The simplest 
cases might be within topological string theory, whose resurgent analysis has also been steadily 
developed [40, 16, 13, 25, 41, 42] in recent years. In particular, the resurgence techniques intro¬ 
duced in [41] have allowed for the generation of large amounts of resurgence data concerning the 
local geometry [42], and are easily extendable to generate equal amounts of data for other 
toric Calabi-Yau geometries such as local P^ x P^. In this context, a very interesting comparison 
is actually now possible due to a recent proposal for obtaining nonperturbative results for topo¬ 
logical strings in toric Calabi-Yau geometries [43] (following up on earlier work in [44, 45]). This 
proposal allows for the calculation of topological string free energies at continuous values of N 
for several local Calabi-Yau geometries; in particular it addresses the examples of local P^ and 
local P^ X P^. It would be extremely interesting to use the aforementioned resurgence data for 
these geometries [41, 42], combined with the methods in the present paper, to investigate how 
accurate would be the match between both approaches, this time around at continuous N. 


2 Exact Finite N Results 


In order to check both the validity and the resulting accuracy of any large N transseries results, 
we need data to compare them against. As already mentioned in section 1, this will be done 
within the gauge-theoretic context of the quartic matrix model, where orthogonal polynomial 
techniques [30] allow for straightforward calculations at small (integer) values of N. We refer the 
reader to [24] for an excellent overview of these methods. 

Let us begin by recalling the definition of the partition function for a one-cut matrix model, 
with polynomial potential V(z), already written in diagonal gauge with eigenvalues Zi (as usual, 
the partition function was originally normalized against the volume of the gauge group U(A^)). 
One has: 


Z{N) 


1 

M 


N 


i=l 




( 2 . 1 ) 
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where A^(z) = Wi^j{zi — Zj)'^ is the Vandermonde determinant. In this expression we have 
chosen the standard convention (within the matrix model literature) of using the so-called string 
coupling, Qs, as the overall normalization, but this is of course just a notation choice and we could 
have likewise used a gauge-theoretic coupling Qs ~ ^ym- case, in the large N ’t Hooft 

limit one has t = ggN fixed and we shall thus mostly work with the ’t Hooft coupling, t, and 
gauge group rank, N, in the following. Furthermore, we shall focus exclusively upon the quartic 
potential with quartic coupling^ A, 

^( 2 ) = ( 2 . 4 ) 

It is a straightforward task to introduce orthogonal polynomials Pn{z) = z"^ + ■ ■ ■, satisfying 

/ dfj.{z) Pn{z) Pm{z) = hn dnm (2.5) 

Jr 


with respect to the measure dp{z) = ff e inherited from the matrix potential. A simple 

calculation [30, 24] then shows how the partition function may be written as a product of the 
above polynomial norms 

N-l N 

Z{N) = n n (2.6) 

n=0 n=l 

or, equivalently, as a product of the recursion coefficients rn = ■ . These coefficients are pre- 

cisely the same as the ones appearing in the recursion relation explicitly constructing orthogonal 
polynomials (for an even potential as we have in our case) 

Pn+l{z) = ZPn{z) - rnPn-l{z). (2.7) 


Conversely, these recursion coefficients may also be written in terms of the partition function: 


rN 


Z{N + l)Z{N -1) 
Z(V)2 


( 2 . 8 ) 


Of course at the end of the day the partition function (2.1) is solely determined by the potential 
(2.4), which means that these coefficients must also be determined by this choice of potential. In 
fact they are; they obey the so-called string equation [30] , a hnite-difference recursive equation 
which in the case of the quartic potential is 


1 - - (r„_i +rn + Tn+i) ) = ngs- 


(2.9) 


®It is important to notice that there are in fact less independent couplings in the problem than apparent at 
first sight. Changing variables as w = leads to 




N\ 


2tv 


( 2 . 2 ) 


with a modified potential 

W (w) = , (2.3) 

where we have defined k = —Xqs. It is then clear that there is a unique independent coupling in the matrix 
integral, although, as we said above, we shall mostly work with the’t Hooft coupling t. 
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Figure 1: Illustration of the (double) monodromy of the confluent hypergeometric function of 
the second kind, for fixed values of a = | and b = ^, and different values of \z\ = 0.1 (lower 
right), \z\ = 1.1 (upper right), and \z\ = 10.1 (left). The figures plot U {a,b \ z) over the complex 
z-plane as argz G (—7r,37r), with the solid line representing the first turn, aigz G (—7r,7r), and 
the dotted line the second, arg^; G (7r,37r). Note the difference in scales between principal and 
secondary sheets, increasingly significant as \z\ grows. In the left plot we have inclosed a zoom-in 
close to the origin, in order to show the trajectory along the first sheet in more detail. 


Note that their relevance extends beyond the above formulae, as these are also the natural objects 
to address and compute when considering the’t Hooft large N limit, see [30]. Furthermore, these 
are also the natural quantities to work with when constructing the large N resurgent transseries 
expansion for the partition function (or for the free energy), e.g., [19, 31]. 

Using these relatively standard techniques, one may now compute the exact partition func¬ 
tion for a few values of the rank N. If N is small, one may proceed and compute the polynomial 
norms, hn, directly. First introduce the moments (under appropriate convergence conditions) 


mn 


dfl{z) z^ 





( 2 . 10 ) 


in terms of which all results that follow may be expressible. 


Of course odd moments vanish as 
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the potential we are considering is an even function. The result is expressible in terms of the 
confluent hypergeometric function of the second kind U{a,b\z) (see, e.g., [26]), which has a 
branch cut along the negative real axis in z, arg 2 ; = ivr, and integral representation 

1 r +00 

U (a,b \ z) = / dx (1 + e“^^, (2-11) 

r(a) Jo 


when Men > 0 and |arg 2 ;| < 
odromy, given by [26] 

[/(a,6|e2"™2) 


Because of the branch cut, this function has non-trivial mon- 


U{a,b\ z) + 


27rie 


—mbm 


sin (Trbm) 
r (1 + a — 6) r (6) sin (vrft) 


M{a,b\ z), m G Z, 


( 2 . 12 ) 

where now M (a,b\ z) is the confluent hypergeometric function of the first kind; an entire function 
with integral representation 


M{a,b\z) 


m 


T{a)T{b-a) 


dx x' 


a—1 


(1 -x)''-“-^e"^ 


(2.13) 


when Me 6 > Me a > 0. In our quartic-potential example (2.10) one always has b = ^, implying 
from (2.12) that we have to rotate twice around the origin in the complex z-plane in order to 
return to the starting point, as illustrated in figure 1. This will also be a distinguishable feature 
for both recursion coefficients, free energy and partition function of our quartic gauge theory. 

One may now explicitly construct the orthogonal polynomials of the quartic matrix model 
in terms of the moments, m^, and find the following first few recursion coefficients 


m 2 

n = —, 

mo 

m4 m 2 

r2 = -, 

m 2 mo 

mo m| - mo m2 
m 2 m^ — m4 mo ’ 

mg m2 m 4 m2 

2 2 
trie tn2 — m| m4 mo — 


mo m2 (mo mo - m 4 m2) 

(mo m 2 - m|) (m4 mo - m|)' 


(2.14) 

(2.15) 

(2.16) 
(2.17) 


Next, it is simple to use these results to find the partition functions at different values of the 
rank. Notice that while the polynomial norms and the recursion coefficients are, in general, 
rational functions, we can see from their structure that certain cancellations occur, ensuring 
that at the end of the day the partition functions are just a sum of products of hypergeometric 
functions. This could have been predicted from the start, since the Vandermonde determinant 
is a polynomial. Explicitly, the partition functions for N = 2, 3,4, 5 are given by 


Z{2) = m 2 mo, 

Z{3) = m 2 (^m4 mo - m^^, 

Z(4) = (me m 2 - m4^ (m^ mo - m^), 

Z(5) = (^m6m2 - m4^ jmg (^m4mo - m^) 


+ mo me mo - m4 m 2 - 


(2.18) 

(2.19) 

( 2 . 20 ) 
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N = 3, I t I = 2, arg (t) e (-TT , 37T) 



Figure 2: Monodromy of the recursion coefficient rs, for different values of |t| = 0.7 (left), 
\t\ = 1.2 (upper right), and |t| = 2 (lower right). The solid line corresponds to the first sheet, 
arg(f) G (—TT,tt), and the dotted line to the second, arg(t) G (tTjStt). In the left plot we inclose 
a zoom-in close to the origin, to see the trajectory along the first sheet in more detail. 


—m4 


me m2 - m4 


)}■ 


( 2 . 21 ) 


Note that for different values of N the above partition functions are not simply numbers, but 
functions. In these formulae they are implicitly functions of the couplings in the matrix integral 
via (2.10), which we trade for the’t Hooft coupling t as usual. As a consequence of (2.12), and as 
we shall now illustrate, these functions display intricate monodromy structures. For instance, in 
figure 2 we illustrate the monodromy properties of rs, for three different values of |t| and varying 
arg(f). Furthermore, the free energy® also displays distinctive monodromy features. The one 
novelty is that due to the logarithm the scale in the second sheet may be considerably larger 

®A word on notation and normalizations: we shall follow standard practice of normalizing the partition function 
with respect to the Gaussian weight (1.5), which means in practice the partition function is in fact given by 


2 : = 




( 2 . 22 ) 


However, with a slight abuse of notation and following common practice in earlier papers, this will always be 
implicitly assumed so that whenever we refer to Z this is always the normalized result Z. The free energy will 
also always be normalized as = log 2, but in this case we shall keep it explicit in notation. 
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Figure 3: Monodromy of the free energy J^(3), when N = 3, for \t\ = 0.6 (left); and of the free 
energy when N = 2, for different values \t\ = 0.8 (upper right) and \t\ = 0.6 (lower right). 

As usual, the solid line corresponds to arg(t) G (—7r,7r) while the dotted line now corresponds to 
arg(t) G ±(7r,7r + 69), where we set 66 = vr/d (left), 69 = Ttt/S (upper right), and 69 = tt (lower 
right). In some plots we have not included the full dotted curves as they show no more relevant 
features beyond what is displayed (their range gets naturally enlarged by force of the logarithm). 


than the corresponding one in the first sheet. This is illustrated in figure 3 for different values 
of N and t (in fact in these plots we do not show the entire range corresponding to the second 
sheet, in order to keep the pictures small enough for illustration purposes). Note that due to the 
logarithm, now the curves do not close upon themselves. 

It is clear that the monodromy features displayed in figures 2 and 3 are quite nontrivial, 
and they change as we change the value of N. What is important to have in mind is that all 
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these features must be precisely captured by the large N transseries which will be introduced in 
the following section. Furthermore, there is one single transseries describing the quartic matrix 
model, and this one transseries must be able to reproduce all these monodromies for all different 
values of N. As we shall see, this will be remarkably captured by instanton physics. 

3 Finite N from Resurgent Large N 

The construction of the perturbative 1/N expansion has a long history. In the matrix model 
context it starts with a (planar) spectral curve [46], out of which one may recursively construct 
the full perturbative series—an endeavor which started out in [47] and culminated in [48]. But, 
as we shall now discuss, there is much more to gauge theory at finite N beyond its perturbative 
1 /N expansion. Depending on the values of the parameters, and on Stokes phenomenon, instan- 
tons may be crucial to achieve exponential accuracy in some results, or, instead, they become 
exponentially enhanced rather than suppressed to completely change the perturbative results and 
correctly reproduce many of the intricate monodromy features we have discussed in the previous 
section. So, how are instantons incorporated into the large N expansion? 

As already explained in the introduction, this is achieved via transseries which go beyond 
the large N expansion by including all its nonperturbative corrections. In practice one deals 
with a formal expression, such as (1.4). But it is important to notice that such transseries 
contain as much information as a would-be analytic expression for whatever function we are 
trying to describe. The explicit connection between the two is achieved by the resummation 
of the former into the latter. This is actually similar to the role that Taylor power-series play 
in describing or representing entire functions. In our case, while similar in spirit, the “power- 
series game” becomes a little bit harder in practice due to the existence of many singularities 
and branch cuts. The prominent role of transseries thus comes about, since finding analytic 
solutions is doomed to fail for most problems, while transseries representations yield natural 
completions to naive perturbative approaches. These in fact include everything needed for a 
complete nonperturbative description of the solutions. Furthermore, treating parameters such as 
the rank N as formal variables removes any previous constraints on their domain of validity {e.g., 
N must be an integer), so that after resummation they can even take real or complex values. 

The observables we shall focus upon are the partition function and the free energy. Let us 
start with the free energy, J^{N, t). In our example of the quartic matrix model, the construction 
of a large N transseries representation for the free energy starts by addressing the coefficients 
rjv (and we shall also display results for these). The large N transseries for rjv is a resurgent 
function of the form 


+ O0 +00 

n{N,t) = (3.1) 

n=0 g=0 

which solves the continuous version of the string equation (2.9) when expressed in terms of the 
’t Hooft coupling. We refer to appendix A, and references therein, for a more detailed explana¬ 
tion of this transseries solution, e.g., the choice of variables, the explicit instanton action, the 
coefficients I3n, the maximum orders up to which we have computed the gf-coefficients, and a few 
explicit examples. Let us nonetheless stress a few points: there are two distinct sums, one in 
the instanton number, n, and the other in the perturbative order, g; the rank N is treated as 
a (continuous) formal parameter; at large N the nonperturbative exponential contributions are 
controlled by the ratio A{t)/t involving the matrix model instanton action (A. 8); the transseries 
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parameter a is so far arbitrary and needs to be fixed for any numerical evaluations; the per¬ 
turbative coefficients grow factorially fast with g, turning each instanton series in 1/N 

asymptotic; and that all building blocks of the transseries are in fact functions of the ’t Hooft 
coupling t. Having determined the transseries for the coefficients r^, one then uses equation 
(A. 12) in order to obtain the transseries for the free energy, 

+ 00 +O0 

T{N,t) = ^9+^-J-W(t). (3.2) 

n=0 g=0 

We once more refer to the appendix for more technical details on this transseries construction. 
Note that when addressing any other gauge theory, the starting point and the method that now 
follows will be exactly the same. The only difference would be a distinct instanton action and 
multi-instanton perturbative coefficients, possibly computed diagramatically, i.e., the difference 
would essentially amount to having a distinct content in the appendix of the paper. 

One thing to notice is that expressions (3.1) or (3.2) are not the most general transseries 
solutions to the quartic matrix model; see [19]. Instead, it turns out that the instanton action 
has a symmetric companion of opposite sign which also solves the relevant differential equation. 
Consequently, the complete transseries solution will depend upon two parameters and it will 
include logarithmic monomials in N associated to resonant sectors. While knowledge of this 
complete transseries is necessary in order to understand the resurgence properties of the free 
energy—namely, how the coefficients of one sector grow and relate to coefficients from other 
sectors—it will not be needed for the particular resummations we address in the present paper. 

Finally, one still has to hx a in order to obtain any numbers at the end of the day. Typically, 
a transseries with a given fixed choice of a will be valid in a specific wedge of the coupling- 
constant complex-plane; in this case of the complex plane associated to the’t Hooft parameter 
t. Different wedges are separated by singular directions on the Borel plane, known as Stokes 
lines. This is where the multi-instanton singularities lie. For our present problem, these singular 
directions on the Borel plane are located at either 0 = 0 or 0 = vr [19, 18], but for the problems 
we shall address in the following we only need to consider the 9 = 0 case. Crossing this Stokes 
line implies that the transseries parameter will “jump” or “turn on”, in the sense that any 
exponentially suppressed contributions previously neglected (as they were invisible behind the 
perturbative expansion) must now be taken into account as they will start growing and eventually 
may take dominance. To make this concrete, consider the free energy transseries (3.2) and write 
it explicitly as J-{N,t,a). Stokes phenomena then translates to the statement^ 

So+R'{N, t, a) = Sq-R'{N, t,a + Si), (3.3) 

where the resummation Sg was defined in (1.3), and where the shift in the transseries parameter 
is controlled by the Stokes constant ^i. A particular case of (3.3) is when we start off with cr = 0, 
meaning the transseries has just the perturbative component, and then after crossing the Stokes 
line the nonperturbative contributions are turned on with a parameter a = Si. 

Now, as we want to bring together and compare the two representations of, say, the free 
energy, the one obtained analytically and the one obtained via resurgent transseries—in fact to 
show that they are equal at integer N —we need to resum the formal transseries into a function. 
Because the transseries is a double sum we must undergo a two-step process which bears the 
name of Borel-Ecalle resummation. The first name deals with each of the asymptotic series 

^In the language of resurgence this jump is captured by the so-called Stokes automorphism; see, e.g., [19, 18]. 
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in 1/A^, while the second takes care of the sum over multi-instantons. Both are important, in 
solving different problems, but, in practice, when we need to get a number out of a particular 
example the Borel resummation takes most of the attention. Let us fix an instanton sector, n, 
and consider the asymptotic series 


+ 0O 

~ J-W(t). (3.4) 

g=0 

As mentioned in the introduction for the perturbative case, Borel resummation first computes 
the Borel transform of the asymptotic series {N, t) —a convergent series in s which may be 

analytically continued—and then evaluates its Laplace transform, yielding the resummation 

/■e®®oo 

SeF^'"HN,t)= I dsB[J^(")](s,t)e"®^. (3.5) 

The choice of the angle 9 for the integration contour must be made carefully due to the sin¬ 
gularities of the Borel transform. Once this is done and all perturbative and multi-instanton 
asymptotic series have been dealt with, one may take the second step and address the sum 
over multi-instantons. This is immediate, so the Borel-Ecalle resummation of the transseries is, 
finally, 

+00 

SF{N,t) = SgF^^\N,t). (3.6) 

n=0 

Note that the left-hand-side is nontrivially independent of 6] all one now has to take into account 
are the Stokes jumps (3.3), i.e., keep track of which wedge in the complex plane are we on. 

When we turn to implement equations (3.5) and (3.6) in an explicit example, such as ours 
for the quartic free energy, it is often the case that we cannot perform the Borel transform 
analytically. This is simply because in most problems nonlinearity prevents us from obtaining 
closed-form expression for the asymptotic coefficients and only a finite number of such coefficients 
are available for computation. The standard approach to circumventing this problem is found in 
using Bade approximants to mimic the analytic continuation of the Borel transform (see, e.g.^ [49, 
26]). Because the Bade approximant is a rational function of s we can capture some of the Borel 
singular behavior. Thus, the Borel-Bade resummation provides a numerical implementation of 
the exact Borel resummation. As such we will define the Borel-Pade-Ecalle resummation of a 
transseries, up to the i-th instanton sector, as® 

SeF^^\N,t) = / dsBB4J^(")](s,t)e-"^ 

Jo 

In the last definition BB^)^^] denotes a (diagonal) order-£ Bade approximant of the Borel trans¬ 
form, and the numerical integration has a cut-off at Smax- The “Ecalle step” of the resummation 
is also truncated in practice, since only a few instanton terms are computed. As we shall see, 
this will nonetheless be more than enough to show the capabilities of transseries resummation. 

®To summarize, we are using S and BP^ to define the numerical approximations to S and B, respectively. 


(3.7) 

(3.8) 
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Figure 4: Approximate complex Borel s-planes for the perturbative (left, and one- 

instanton (right, free energies obtained by plotting poles of the Fade approximant when 

t = 2. Due to limited computational resources we have less points for higher instanton sectors, 
as compared to the pertnrbative sector. Still, the instanton action singularities are very clear, 
with the accumulation of poles signaling their associated logarithmic branch cuts. 
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Note that what we have described is not a numerical method, but rather a numerical approx¬ 
imation to an analytical procedure. As we already mentioned, the Borel-Ecalle representation 
of a function is somewhat analogous to the Taylor power-series representation of (another) func¬ 
tion, and this is what (3.5) and (3.6) set up, out of the transseries. Of course if one is to extract 
a number ont of any of these analytical representations, some approximation (or truncation) 
mnst be considered. In Taylor power-series one jnst trnncates at a given order and then sums. 
In the present Borel-Ecalle framework, where the functional complexity is larger, one needs to 
implement the above Borel-Pade-Ecalle resnmmation. 

We can see an example in figure 4. The singularities on the Borel-Pade plane convey the 
image of branch cuts, where the branch points are given by the instanton actions that appear 
in the transseries. As we have commented before, instanton actions come in pairs of opposite 
signs. In here, we also notice the presence of a displacement of the instanton action A{t)/t by 
a constant term 27ri. This is in agreement with general expectations of [50], where instanton 
actions in matrix models and topological strings should be linear combinations of spectral curve 
B-periods, such as in [16, 17], with spectral curve A-periods, such as in [25] (i.e., the factor of 
27rit). Note however that due to the nature of the string equation, the sector associated with 
this other action is indistinguishable from that of A(t)/t alone, and the two can be combined. 
In fact if we were to consider a two-parameter transseries® with sectors associated to A and 
A = A± 27rit, denoted by [n\n], then after plngging such an ansatz into the string equation we 
would find that the “mixed” coefficients satisfy 


r>ln\n] 

g 


(n + n) 
n! n! 




(3.9) 


Plngging this resnlt back into its two-parameter transseries wonld reduce it to the one-parameter 
transseries (3.1), with a simple shift in parameter a —)• (1 -|- e^^^'^)(T. In the case of integer N 

®This should not be confused with the two-parameter transseries in [19], where the two actions are ±A. The 
inclusion of this second action A = Ai: 2nit in that set-up would lead to a four-parameter transseries. 
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we have considered so far, the factor is equal to 2. But since the transseries parameter a needs 
to be fixed in any case (which we will do next), the effect of the second instanton action with 
the 27ri displacement is already automatically included in our results. 

The final issue to address concerns the fact that equations (3.1) and (3.2) are actually 
representing a family of transseries, indexed by the transseries parameter a. The resummation 
procedure cannot be complete until a is specified. As we have discussed earlier, its particular 
value is subject to Stokes transitions that may add or subtract the Stokes constant Si, but it 
still needs to be fixed at some point, with Stokes transitions then specifying it wherever else. In 
the present section we are focusing upon the case t G M'*', where one can numerically check that 
the transseries parameter is 

(3.10) 

equal to the Stokes constant Si [16, 13, 19]. Let us recall that A is the quartic coupling-constant 
which we are going to set to —1 without loss of generality. More complicated gauge theories in 
more “physical” scenarios may eventually require that the fixing of transseries parameters must 
be done against some laboratory measurement. 

Now that all the ingredients are on the table, we can explicitly show how adding more and 
more instanton contributions of the transseries gets us closer and closer to the exact result. In 
table 1 we display explicit numerical examples for the recursion coefficients (left) and the free 
energy (right). The first four rows correspond to the BoreBPade resummations (3.8) of each 
sector, the fifth row is their sum (3.7), and the sixth has the analytical results, (2.16) and the 
logarithm of (2.19). In the total result we have labelled with different colors^*^ the digits of 
the exact result that are matched after including the perturbative (blue), 1-instanton (green), 
2-instanton (yellow) and 3-instanton (red) sectors. One can clearly see how the instanton contri¬ 
butions^^ provide exponentially small corrections with respect to the perturbative contribution, 
where the size of the correction is naturally controlled by the instanton action. 


Sector 


Perturbative 

2.615 

796 

570 

569 

705 

50... 

-1.973 

899 

279493 

161 

74... 

1-Instanton 

0.000 

487953 

495 

567 

22... 

-0.000 

020 

359 

080 

917 

15... 

2-Instanton 

0.000 

000 

009 

807 

788 

15... 

-0.000 

000 

000 

300 

789 

88... 

3-Instanton 

0.000 

000 

000 

000 

245 

38... 

-0.000 

000 

000 

000 

004 

71... 

Total 

2.616 

284 

533 

873 

306 

27... 

-1.973 

919 

638 

874 

873 

50... 

Exact 

2.616 

284 

533 

873 

306 

26... 

-1.973 

919 

638 

874 

873 

50... 



Table 1: Comparison of the truncated TZ and J- transseries, up to the instanton sector n = 
0,1,2,3, against the exact result for t = 6 and N = 3. All digits displayed are stable. On the 
left table the last digit of the resummation must be corrected by the 4-instanton contribution. 


We illustrate this visually in figure 5, also for different values of the rank. On the x-axis we 


vary N and on the y-axis we plot — log;^o 


- rN 


(left) and — log;^o — T'(A^) (right). 


'^°This is in fact the color code we shall be using in all subsequent plots. 

'^^There is a slight abuse of notation in the column titles when we write them as and What we 

show is actually the whole contribution to the resummed objects, as in (3.7), so that herein SqTZ^''' and are 

in fact multiplied by the relevant powers of cr and exp{—NA{t)/t). 
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N N 


Figure 5: Number of decimal places up to which the resummation of the recursion coefficients 
(left) and the free energy (right) match their exact counterparts, with N = 1,..., 5 and t = 6. 



N 


Figure 6: Resummation of the partition function for continuous N G (0, 6), t = 1, including up 
to 3-instantons (the lines are all superposed), and the exact results at integer values of N. 


where rjsr are the exact recursion coefficients (2.14)-(2.17) and iF{N) are the logarithms of the 
exact (normalized) partition functions (2.18)-(2.21). This quantity effectively tells us the number 
of decimal places to which the analytical result and the resummation agree, and as we saw in the 
example above we get closer and closer to the full answer as we add more and more instanton 
sectors. One can check that at = 3 the matched digits are the ones shown in table 1. Note 
that we have chosen a relatively high value of t where the picture is clearer. As we move to lower 
t the instanton contributions get smaller in absolute value, and it may happen that if has 
not stabilized at enough digits it becomes harder to see their effect. 

Having seen how the transseries so precisely captures the exact results, at small integer 
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values of the rank, one may ask if it can go beyond this requirement and actually compute the 
free energy or partition function at continuous values of A^. It should be clear that nothing 
changes as one considers the resummation for non-integer N. Furthermore, as we have analytical 
results for the exact partition function at integer values of N, we can ask how the resummation 
interpolates between them. This is shown in figure 6, with t = 1, where it is clear that the 
resummation produces a smooth interpolation between the analytical results arising from the 
matrix integral. Of course these results are already going well beyond the exact matrix integral 
results, as the latter are not even defined when N is non-integer. 

It is worth pointing out that we are not resumming Z in the same way we did for iF or TZ, 
since it is very inefficient to exponentiate the free energy transseries and then extract the 1/N 
coefficients at each order. Instead, we first resum the free energy and only then exponentiate the 
result. Furthermore, note that for the almost entirety of the plot in figure 6, the resummations 
including any or all instanton sectors are indistinguishable and the four lines are coincident. 
As we showed above, these distinctions only appear after a certain number of decimal places, 
and this is impossible to spot in these scales. However, in the next section we will also explore 
complex-valued t and we will see cases where there is a “macroscopic” distance between different 
resummations. Finally, we notice that the coincident lines seem to split apart as we get closer 
to = 0. The fact that we are dealing with a normalized partition function means that 
Z{0) = 1. However, this is a point of infinitely strong coupling, and even the large amount of 
data we have is insufficient, from a numerical standpoint, to get consistent results at this point 
(i.e., technically, at infinitely strong coupling we would need infinite terms in the Borel-Ecalle 
transseries resummation, which is unachievable). 

4 Analytic Continuations and Stokes Phenomenon 

In the previous section we have limited our attention to the case where t G ]R"*“. While definitely 
crediting the power of resurgent transseries in achieving to go beyond integer rank and actually 
define the gauge theory partition function at continuous values of N, the reader might get 
the wrong impression that other than that all the transseries is implementing is smaller and 
smaller instanton corrections to the perturbative expansion, piling on top of each other. This 
is certainly not the generic case and such picture will dramatically change as we consider the 
analytic continuation onto complex values of the’t Hooft parameter, due to Stokes phenomenon. 

To understand how this works, let us first address an illustrative example. Consider the 
free energy when N = 3 and at fixed’t Hooft coupling |t| = 0.6 but varying argument arg(t) G 
(—7r,7r -|- 7r/4). This is a continuous function with nontrivial monodromy, as shown in figure 3, 
computed from analytic expressions. In the ’t Hooft limit, the large N resurgent transseries 
for iF, (3.2), must have all the required information in order to reproduce this plot. To extract 
it we have to resum it, as explained in the previous section, but also take into account Stokes 
phenomenon. This last step is crucial and in practice it implies selecting a particular member 
from the family of transseries parametrized by a. For t G M"*" we saw that the correct value is 
given by (3.10), but as we now vary arg(t) we will need to implement Stokes transitions as shown 
in (3.3). This is a general feature of representing functions as transseries: we need the whole 
family of solutions, as parametrized by the parameter a, and a practical understanding of how 
Stokes phenomenon selects the right member as we move across the complex t-plane. 

Stokes transitions take place at Stokes lines, the rays where multi-instanton singularities lie 
on the complex Borel plane. In the complex ’t Hooft plane they satisfy Im(A(f)/t) = 0, and 
impose the familiar jump in the transseries parameter. One may also cross anti-Stokes lines, 
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Figure 7: Phase diagrams showing anti-Stokes (thick, dashed, blue curved) and Stokes (thick, 
solid, red curved) boundaries along with the path |t| = 0.6 (arrowed black) and the intersections 
of this path with Stokes (thin, solid, red) and anti-Stokes (thin, dashed, blue) lines. The left 
plot represents arg(t) G (—vr,-|-7r) while the one on the right is for arg(t) G (tt, Svr). In the last 
line of the caption we show how the motion on the two sheets takes place as arg(t) is changed 
from 0 (the solid square on the first diagram) all the way (back) to dvr. Note that there is a third 
intersection with an anti-Stokes line in the second sheet, for arg(t) = 27r, but we do not draw the 
intersecting ray because it would lie along the positive real axis. 


where Me{A{t)/t) = 0, with instantons taking dominance over the perturbative series—in fact 
at this point all contributions, perturbative and nonperturbative, are of the same order. As we 
move beyond it and into a region where Me {A[t)/t) < 0, the instanton contributions are no longer 
exponentially suppressed. However, this does not mean that the transseries representation breaks 
down. The transseries is a formal object that includes complex instantons in general, but where 
the label “exponentially suppressed” or “exponentially enhanced” only applies in the formal large 
N limit. In the resummation process where N becomes finite, even small, and where we may 
venture into the complex plane, this distinction is somewhat irrelevant and even not appropriate 
any more. The separation between perturbative and (multi) instanton sectors is only set up at the 
initial definition of the transseries, where one finds out that for some values of t the resummation 
of the perturbative series alone is enough to give the full answer, thus having cr = 0. But this 
washes aways at finite N and complex t where the resummation of the transseries yields, in 
practice, a power-series in a. All these features will be clear in the examples that follow. 

Graphically we can represent the Stokes and anti-Stokes boundaries in the complex t-plane^^ 
and then, for each particular case, determine the corresponding lines as intersections with the 
path \t\ = constant. We show this in figure 7 for our example value |t| = 0.6, where the diagram 
on the left corresponds to the first sheet, arg(t) G (—vr,vr), and the one on the right to the 
second sheet, arg(t) G (tt, Svr). We see that in the first sheet there are four Stokes lines at 

'^^These become the large N phase diagram for the quartic matrix model; see [51-53] for further discussions. 
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N = 3, I t I = 0.6, arg (t) e (-7T , 57T/4) 



-•— Perturbative 
-•— + 1-instanton 

-•— + 2-instantons 

-•— + 3-instantons 

S — Exact 

- Stokes Lines 

- Anti -Stokes Line 

- An. Cont . Border 


Figure 8: Resummation of the free energy transseries for |t| = 0.6 and N = 3, taking into account 
the different values of a (due to Stokes phenomenon) depending on arg(t) G (—tt, vr + vr/d). The 
Stokes lines illustrate how in some regions a = 0 and the perturbative resummation yields the 
correct result, while in other regions a = Si and we need to include up to three instantons in 
order to reproduce the exact result (shown in the left plot of figure 3 and enclosed here). It is 
also clear how past the anti-Stokes phase boundary the perturbative contribution is no longer 
reliable and instantons are required to take over in order to yield the correct monodromy results. 
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N = 3, I t I = 0.7, arg (t) s (-rr , 57T/4) 


- Planar 

-•— Perturbative 

• + 1-instanton 

— + 2-instantons 

-•— + 3-instantons 

■©— Exact 

- Stokes Lines 

- Anti -Stokes Line 

- An. Cont . Border 


-4 -3 -2 -1 0 

Re ^{3) 

Figure 9: Resummation of ra, for \t\ = 0.7, in direct correspondence with the left plot in figure 
2. Note the various Stokes lines in the principal domain for arg(t), and the anti-Stokes phase 
boundary in the second domain indicating that instanton corrections have become of order one. 



different arguments arg(t) G (—7r,7r), and in the second sheet there are three anti-Stokes phase 
boundaries. There is also a symmetry with respect to the real line, so we can focus on arg(t) > 0. 
Of course the information shown in figure 7 is incomplete without the actual resummation of the 
transseries, for the different values of arg(t), and that is displayed in figure 8. In the following we 
shall explain both these plots^^ at the same time, moving along the arrowed path drawn on the 
phase diagram, figure 7, and then looking at the relevant features of the resummed transseries, 
figure 8. 

For arg(t) = 0 we are in a case similar to the one described in the previous section: we can 
reproduce the value of the free energy with increasing accuracy by piling up smaller and smaller 
instanton contributions on top of the (already quite accurate) perturbative result. The value of 
the transseries parameter is (3.10), i.e., the Stokes constant. The reader may also want to take 
another look at table 1 and figure 5 from the previous section^"^. As we increase arg(t) from 0 
onwards we will cross the first Stokes line, where a transition occurs that selects the transseries 
with cr = 0. That is, we reach a region where perturbation theory alone is enough to reproduce 
the exact value of the free energy. To show that this is the case, we display the different instanton 
contributions in table 2. It is clear that if we were to add them to the perturbative result we 
would immediately deviate from the exact result! Pushing arg(t) further towards +71 we cross the 

very pedagogical introduction to this type of Argand plots may be found in the excellent review [54], 
addressing the simple example of the Airy function (which is a linear problem with no multi-instanton sectors). A 
similar example but for the Bessel function (again linear and without multi-instantons) appears in [55]. 

^^Note that in the present section we focus on a smaller value of ]t] than in the previous section, to illustrate 
Stokes phenomena. However the precision of the results is reduced, due to the finite number of available transseries 
coefficients, fn practice, the contributions of the second and third instantons do not provide stable digits to display. 
In any case, in figure 8 the different contributions cannot be distinguished with the naked eye. 
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N = 3, I t I = 1.2, arg (t) s (-157T/8 , 157T/8) 


- Planar 

-•— Perturbative 

-•— + 1 -instanton 

-•— + 2-instantons 

-•— + 3-instantons 

O Exact 

- An. Cont . Border 


-4 -3 -2 -1 0 

Re 9?(3) 

Figure 10: Resummation of r^, for |t| = 1.2, in direct correspondence with the top right plot 
in figure 2. This plot illustrates very clearly how the perturbative resummation keeps following 
the planar approximation (and thus producing an incorrect result), while instanton effects take 
dominance in order to produce the correct monodromy results. 


Sector 

Me 


Perturbative 

+0.130 991945 237228... 

n nnn n^n a^a ^KnaAA 

-0.478 840 360187836... 

n nnn ni n c^^n nov 

i Xllo 

n nnn nnn nno nn7 

KJ ,\JKJKJ \J J.\J (JKjyj tJLJ i «J\JO... 

n nnn nnn nni 97q 907 


n nnn nnn nnn nnn nnv 

\j.yj\jyj uuw uux o t o i . . . 

n nnn nnn nnn nnn nnc: 

llio LiCtlJ-L vjli 

Exact 

+0.130 991945 237 228 ... 

-0.478 840 360187836... 



Table 2: Comparison of the real and imaginary parts of the resummed J^-transseries, at the 
instanton sector n = 0,1, 2, 3, and compared against the exact result for N = 3 and t = | 

All digits displayed are stable. Note that all the digits in the perturbative resummation already 
match the exact solution, so the transseries parameter must be 0. 


second Stokes line. This restores the value of a back to what it was at arg(t) = 0, namely (3.10). 
After this point the perturbative result alone is already significantly different from the exact one, 
as can be clearly seen in figure 8. We could, in principle, stop at arg(t) = +7r, but we know that 
the partition function has monodromy 2 around t = 0 so we keep pushing the calculation. In 
order to do this, for the exact value of the free energy we will use the analytic continuation of the 
hypergeometric function, explained around equation (2.12). For the transseries we need to do 
an analytical continuation of both its coefficients and the instanton action. As we move beyond 
arg(t) > +7r the first observation is that perturbation theory becomes less and less accurate. 
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N = 3, I t I = 1.2, arg (t) € (-7T , 3 tt ) 



Planar 
Perturbative 
+ 1-instanton 
+ 2-instantons 
+ 3-instantons 
Exact 

An. Cont . Border 


Figure 11: Resummation of rs, for \t\ = 2, in direct correspondence with the bottom right plot 
in figure 2. Note how the perturbative resummation keeps following the planar approximation. 


N = 2, I t I = 0.8, arg (t) e (-7T , 77 t/4) 


• Perturbative 

• + 1-instanton 
♦- + 2-instantons 

• + 3-instantons 
O Exact 

- Stokes Lines 

. Anti -Stokes Line 

. An. Cont . Border 


0 2 4 6 8 

Re r{2) 

Figure 12: Resummation of for \t\ = 0.8, in direct correspondence with the top right plot 
in figure 3. While the perturbative resummation still follows the general trend of the monodromy, 
it is no longer reliable past the anti-Stokes phase boundary. 



Shortly afterwards we cross an anti-Stokes line, signaling a change in the sign of Me {A(t)/t), from 
negative to positive. Now instanton contributions are roughly of the same order of magnitude as 
the perturbative contribution, and this effect is dramatically enhanced as arg(f) grows. Following 
the sequence of points in figure 8 we see how each new instanton contribution struggles to move 
the resummation line closer to the exact curve. As we get to arg(f) = -|-7r -|- ^ not even three 
instanton terms can give an accurate result, and higher instanton sectors are needed to keep up 
with the analytical curve. 


- 22 - 



















Re !r(2) 


• Perturbative 

• + 1-instanton 

— + 2-instantons 

• + 3-instantons 

■©— Exact 

. An. Cont . Border 


Figure 13: Resummation of -F(2), for \t\ = 2, in direct correspondence with the bottom right 
plot in figure 3. Unlike previous plots, note how here the perturbative result is pretty reliable. 


Figure 8 is a very rich picture which involves both Stokes and anti-Stokes lines, regions 
with fj = 0 and with a ^ 0, and a crisp image of the importance of instanton corrections. In 
fact it shows how the instantons in the transseries are relevant for much more than achieving 
exponential accuracy in matching the exact results: they actually need to take dominance in 
order to properly describe the physics at small N. This should be extremely compelling evidence 
towards the relevance of resurgent transseries in describing the nonperturbative realm. We may 
now proceed with exploring the gauge-theory parameter space. In figures 9, 10, and 11, we 
expand our list of examples by keeping N = 3 but varying the value of |t|; and this time around 
we consider the orthogonal polynomial recursion coefficients TZ{N,t). In some of these figures 
for TZ{N,t) we can actually follow the entire monodromy and “close” the curves (the cases we 
are presently addressing are the ones shown earlier in figure 2). These figures also show a rather 
interesting feature: the perturbative resummation has a tendency to follow along the planar 
approximation and thus, past anti-Stokes boundaries, it completely misses the correct features of 
the problem. It is the instantons that come to save the day and properly describe the monodromy 
properties we are trying to reproduce. Furthermore, we keep expanding our list of examples in 
figures 12 and 13, by now addressing the free energy at = 2 and for new values of |t|; and 
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N 


Figure 14: Real part of the resummed partition function, SqZ, for continuous negative to 
positive N G (—6, 6) and fixed t = \ e“^. The accuracy is no longer reliable for N < —4. 



4 -0.5 


(3) 

Figure 15: Real part of the partition function, resummed up to three instantons, Sq for 

QOtt 1 

complex rank Me(iV) G (0,4) and Im(A^) G (—1/2,1/2), and fixed t = lOe^w. 


these cases are in correspondence with the exact results shown in figure 3. In all these figures 
we can cleanly identify all the Stokes, anti-Stokes, and analytic continuation transitions which 
occur as we draw the constant |f| paths on the phase diagram of figure 7. 
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Finally, we wish to study the partition function as a function of N, where N will be taken as 
an arbitrary complex number, and at fixed ’t Hooft coupling t. The first thing one can try to do 
is extend our results towards negative N. This is shown in figure 14, for t = 5 e“^. At positive 
N the different lines smoothly interpolate between the exact results, and it is in fact impossible 
to tell them apart with the naked eye. Only very close to = 0 do we notice the instanton 
contributions behaving incorrectly. As we have discussed earlier, this is a point of infinitely strong 
coupling where we do not have enough data—in fact at this point the perturbative resummation 
is more reliable as we have more than twice the perturbative coefficients as compared to the 
(multi) instanton coefficients (see the appendix). As one moves towards A” < 0, we should bear 
in mind that we have entered a region where exp (—NA(t)/t) has changed sign. This means 
that the perturbative sector is no longer a viable approximation to the full answer, since it just 
became of the same order as the instanton sectors. The plot in figure 14 illustrates how up until 
N ~ —4 the instanton sectors still look consistent, and are providing corrections which sit on 
top of the previous sector. However, as N gets increasingly negative we start getting wilder 
oscillations and the resummation can no longer be trusted. We would need more coefficients, 
and crucially higher instanton sectors, in order to carry on. In spite of this, it is clear how the 
transseries resummation allows us to define gauge theory at negative rank! 

90771 

We can next look at complex values of N, and this is shown in hgure 15, for t = lOe^w (in 
this figure we are just plotting the third instanton contribution to the partition function, which 
visually is in fact indistinguishable from the perturbative or lower instanton resummations). Sim¬ 
ilarly to what we have found when looking at real N, positive or negative, we have a function 
which oscillates and those oscillations can get milder or harsher as we move away from the real 
line. Note that in [27] the authors were able to analytically show that the M = 8 ABJ(M) parti¬ 
tion function was an entire function of N. Crucial to that, as we have mentioned before, was the 
drastic simplification they found in their nonperturbative structure. The quartic matrix model 
partition function we address in the present paper has no supersymmetry, milder integrability 
properties, and a full transseries to deal with, implying that such an analytical proof is at present 
still unachievable^^. Generic gauge theories will be even more complicated. Nevertheless, we can 
carry out numerical tests, such as the ones illustrated above, and based on the evidence we have 
been able to produce we find encouraging signs that the partition function will indeed be an entire 
function of N, at hxed ’t Hooft coupling. Future research should definitely address this issue, 
perhaps starting by producing more data, both higher 1/N coefficients in the sectors we have 
already computed and in higher instanton sectors. In any case, it is certainly remarkable that 
our one single transseries was able to correctly reproduce the many different (nonperturbative) 
monodromy structures at play, for different values of N and t. It is also extremely interesting 
that at the same time it went far beyond any available analytical results by extending the gauge 
theory to arbitrary complex N. 
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'^®Further studies of the partition function of the quartic matrix model will appear in [52, 53]. 
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A Quartic Matrix Model: Transseries Data 

Our resurgence approach to extract finite N results out of the (nonperturbative) large N expan¬ 
sion was illustrated throughout the main body of the paper within the example of the transseries 
solution to the one-cut quartic matrix model. In this appendix we briefly overview the data 
for this gauge theory, and the methods required in order to obtain it. These methods were ex¬ 
tensively used in [19, 31] to construct full two-parameter transseries solutions to the one- and 
two-cut quartic matrix model, resulting in large amounts of data for both perturbative and multi- 
instanton sectors. For the purposes of our current numerical analysis, however, the full data in 
[19, 31] is not needed but just about “half”. The transseries written down in [19, 31] were written 
as asymptotic expansions in the matrix-model string coupling, gg, being a bit closer in spirit to 
the string theoretic literature. In this case, remaining closer to the gauge theoretic literature, we 
shall slightly rewrite these results using the’t Hooft coupling instead, where t = QsN. 

The starting point is the so-called string equation, an equation computing the recursion 
coefficients r„, (2.6), given some choice of polynomial potential in the matrix integral for the 
partition function. As discussed earlier, if one is able to compute these coefficients, then the 
partition function itself follows. For the case of the quartic potential, the string equation takes 
the form 

rn(^l-^(rn-i+rn + rn+i)^ =ngs. (A.l) 

The standard procedure starts by considering a continuum limit, where we introduce a new 
continuous variable x = ngs and where the recursion coefficients get promoted to functions 
TZ{x) (we shall take the value x = tin the following). This function is then written as a transseries, 

+00 

TZ{N,t) = ^^T"7^('^)(iV,^), 

+ 00 

9=0 

where a is the transseries parameter and A{t) is the instanton action. The coefficient fin may 
be regarded as a “characteristic exponent”. In the present case fin = n/2 (but see [19] for 
tables indicating the different values of fin for the transseries of the quartic model). Note that 
in general one actually needs to consider a two-parameter transseries, with instanton actions 
izA(t), leading to generalized instantons as one needs to consider many new nonperturbative 
sectors (n) —)> (njm). Moreover, all these nonperturbative sectors are related to each other via 
resurgence; see [19, 31]. For most of the numerical analysis we carry out in this paper, the above 
one-parameter transseries is enough. Another thing to notice is that the instanton action and 
the asymptotic coefficients have “attached” adequate powers of t. This is just so that what we 
herein call A{t) and R^^\t) are the exact same quantities as those computed in [19]. Now, by 
plugging this transseries expansion into the (continuous) string equation^® 

n{t) 11 - ^ (7^(t - 5.) + n{t) + n{t + gg))^ =t, (a.4) 

we can recursively solve for the coefficients (t ), as well as compute the instanton action (see 

[19] for details). The equations one finds are differential when re = 1 and algebraic in all other 

^®Not to clutter notation too much, we sometimes omit the first argument in TZ(N,t). 


(A.2) 

(A.3) 
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cases. The coefficients R^\t) can be written in terms of the variable r = ^ (l — \/l — 2Xtj and 
they have a pattern of the form 




(Ar)Pi 


rP2 (3 - 3A rf'-^ (3 - A ^ 
The exponents in the prefactor are fnnctions of n and g, 


P<”l(r). 


Pi = \ - 2 ), 

P 2 = n + g - 1, 

P3 = \ + lOg - 4), 

P 4 = 1 (3n + 6g + 26-4), 


(A.5) 


(A.6) 


with 6 = n mod 2, and the Pg^\r) are polynomials of degree (n + 6g + 6 — 2)/2. For the pnrposes 
of the numerical analysis in the main text, we have focused only on instanton sectors up to n = 3. 
The maximum order of the asymptotic coefficients computed in each sector is shown below, in 
table 3. For completeness let us also write down the first coefficients in each sector. For n = 0, 


n 

0 

1 

2 

3 

^max 

200 

50 

50 

50 


Table 3: Highest order g for which we have calculated R^\t). 


= r, 


o(0) _ A r 
K 2 — 


6(1- Ar) 


4 ’ 


^( 0 ) _ 7AV (5 + 2Ar) 
— 


72(1 - Xry 


(A.7) 


With n = 1 the equation at order N ^ gives a solution for At order we compute 


the instanton action (see as well [16, 13, 17]) 

/3 — 2Xr\ 1 

A = -- (2 - Ar) arccosh i ——— j + — v^(3 - Ar) (3 - 3Ar), 


and the first coefficients that follow are 


dW _ 

Rq — 


\/W 


(3-Ar)^/^ (3-3Ar)^/^’ 
Finally, for n = 2,3, we have 
A^r 


0 ( 2 )__ 

Rq — 


0 ( 3 ) _ 

Rq — 


2(3-Ar)^/'^ (3-3Ar)^/^’ 
3 (Ar)^/^ (2 - Ar) 


8r2 (3-Ar)^/^ (3-3Ar)^^/^ 


( 1 ) O-y/W (6 — 3Ar — 6A^r^ + 2A^r^) 


Ry> = - 


8r (3 - Ar)^/^ (3 - 3Ar) 


ni/4 


( 2 ) _ 3A2 (18 + 117Ar - 102AV2 + 22X^r^) 


R}{’ = 


8 (3 - Xry (3 - 3Ar)" 


( 3 ) 27 (Ar)^/^ (18 + Ar - 12A2r2 + dA^r^) 

, Ri — 


16r3 (3 - Xrf^/^ (3 - 3Ar)^®/^ 


(A.8) 


(A.9) 


(A.IO) 

(A.ll) 
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Having addressed the transseries multi-instanton structure of TZ{N,t), one may proceed and 
address the free energy next. The relation between these quantities is encapsulated in a Toda-like 
equation (which is in fact nothing more than the continuous version of (2.8)), 

T{t + gs) - 2T{t) +T{t-gs) = log . (A.12) 

Here T = F — Fq denotes the free energy of the quartic matrix model normalized against the 
Gaussian contribution. By force of the above relation, the free energy will inherit the transseries 
structure of TZ{N,t), such that again one finds 


with 


+ 00 



F{N,t) = 


(A.13) 


+00 

-2g^2g-2^g)(^)^ 


t) 

9=0 

(A.14) 


+ 00 


t) 


(A.15) 



9=0 



where = n/2. Plugging this back in (A.12) and expanding in powers of a we find 


J-(-) (t + ffs) - (t) + {t-gs) = (t), 


(A.16) 


where <l>^”)(t) is the n-th instanton sector of log ■ The hrst few sectors, which we will be 

using, are 


chW(t) = log , 


$W(t) 


i?W(t) 

i?(0)(t)’ 


$(2)(t) 




1 / 

2 yRio){t) 


2 

. (A.17) 


The standard course of action, e.g., [30, 13, 19, 31], essentially amounts to inverting (A.12) and 
finding explicit equations for each transseries component (t). Herein we have taken a slightly 

different route, already starting at the perturbative level, n = 0, which turns out to be much 
more computationally efficient at high genus. What we do is simply to use (A.12) directly: start 
by expanding its right-hand-side in powers of 1 /N, where the expansion of the logarithm is now 
written as a sum over partitions 


log 



log 




+ 00 

+ E 


k=l 


]\[2k 


^(-1) 

s>l 


s-1 




E 


<’(*) 

R'ihti 


(A. 18 ) 


Now on the left-hand-side, at order we see that the highest order coefficients cancel 

out and we are left with a differential equation for d^F 2 ^\ Knowing that the perturbative free 
energies follow a pattern^^ written in terms of the variable r (see [19]), 


-^25 


^^With exceptions at g = 0,1. 


A29-1 

( 2 -Ar)23-2(1 _ 



(A.19) 


- 28 - 














with (r) a polynomial of degree 3g — 3, it is a computationally straightforward task to 

plug this ansatz into the equation and solve for the coefficients of the polynomial. The boundary 

conditions amount to setting = 0, which is a consequence of the Gaussian normalization. 

^ A=o 

The highest order to which we computed T”® is shown in table 4. The first coefficients are 


•7r(0) 

-^0 


— Xr (9Ar - 16) + 12 (2 - XrY log 
96 


2 - Ar 


-riO) 

•^2 

-p{0) 

-^4 


-log , 

12 ^ V 2 - Ar 

A3 (4i_x2^3 _ i 85;\^2 200r) 

2880 (Ar - 2f (1 - Xrf 


(A.20) 

(A.21) 

(A.22) 


Proceeding with the one and two instanton coefficients, the calculation follows in a straighforward 
fashion. The equations at order Qs are now algebraic equations for and we solve them 

using the same method of plugging-in an ansatz and reducing the problem to one of finding 
coefficients of a polynomial. The maximum orders of the coefficients we found for n = 1,2,3 are 
shown in table 4. Essentially, we are limited by the number of Rg ^ we calculated beforehand. 
The functions J-g"'\ which we computed for n = 1, 2, 3, are of the form^® 


n 

0 

1 

2 

3 

^max 

130 

50 

50 

50 


Table 4: Highest order g for which we have calculated Rg^\t). 


jr{n) 
'' 9 


(Ar)Pi+i 

rP2+i (3 _ (3 - Ar)^^+^"'^ 



n > 1, 


(A.23) 


where the 'Pg^\r) are polynomials of degree {6g n — S)/2, and the exponents pi were defined 
in (A.6). At lowest order, the n = 1, 2 and 3 coefficients are 


F 

f 


( 1 ) 


•r(l) 

•^1 

•77(2) 

•^0 

-p{2) 

•'I 

-p{3) 

■^0 


A3/2 ^1/2 

2(3-3Ar)'^/^ (3-Ar)^/^’ 

9A3/2 (0 75Ar _ 54A2r2 IOAV^) 

16ri/2 (3 _ 3Ar)i'5/^ (3 - 
_ AV(3-2Ar) 

8 (3 - 3Ar)®/^ (3 - ’ 

3A3 (54 + 531Ar - 846AV2 4Q2X^r^ - 92A^r^) 

32 (3 - 3Ar)^ (3 - Xrf 

A 9 / 2 ^ 3/2 ^-0 _ 

48(3-3Ar)^®/^ (3-Ar)"^/^’ 


(A.24) 

(A.25) 

(A.26) 

(A.27) 

(A.28) 


note that there is a small typo in the formulae for the free energy coefficients in [19]. The factors of t — 
therein, where c? = r, should in fact be o? — t. 
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(A.29) 


( 3 ) _ ('_;^o8 - 1044Ar + IDlTA^r^ - lOSOA^r^ + 212A^r^) 

128 (3 - 3Ar)^®/'^ (3 - 

Finally, let us note that even though we have presented several results for the partition func¬ 
tion in the main body of the paper, we have actually not addressed its transseries representation. 
In fact, it is computationally very inefficient to exponentiate the free energy transseries and then 
extract the transseries coefficients for the partition function. If we had done that, we would 
have ended up only going to much lower orders than what we have achieved for the free energy. 
As such, we have instead always performed any required resummation first, for the free energy 
transseries, and only then exponentiated the result. 
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